Hands-on Exercise 3b - Programming Animated Statistical Graphics with R

3.1 Overview

In this exercise, we will learn how to create animated data visualization by using gganimate and plotly r packages. We will also learn how to reshape data using tidyr package and process, wrangle and transform data using dplyr package.

3.2 Getting started

3.2.1 Loading the R packages

First we install and load the following R packages:

• plotly, R library for plotting interactive statistical graphs.

• gganimate, an ggplot extension for creating animated statistical graphs.

• gifski converts video frames to GIF animations using pngquant’s fancy features for efficient cross-frame palettes and temporal dithering. It produces animated GIFs that use thousands of colors per frame.

• gapminder: An excerpt of the data available at Gapminder.org. We just want to use its country_colors scheme.

• tidyverse, a family of modern R packages specially designed to support data science, analysis and communication task including creating static statistical graphs.

pacman::p_load(readxl, gifski, gapminder,
               plotly, gganimate, tidyverse)

3.2.2 Importing Data

In this hands-on exercise, the Data worksheet from GlobalPopulation Excel workbook will be used.

Importing Data worksheet from GlobalPopulation Excel workbook by using appropriate R package from tidyverse family.

col <- c("Country", "Continent")
globalPop <- read_xls("data/GlobalPopulation.xls",
                      sheet="Data") %>%
  mutate_each_(funs(factor(.)), col) %>%
  mutate(Year = as.integer(Year))

Note

• read_xls() of readxl package is used to import the Excel worksheet.

• mutate_each_() of dplyr package is used to convert all character data type into factor.

• mutate of dplyr package is used to convert data values of Year field into integer.

mutate_each_() was deprecated in dplyr 0.7.0. and funs() was deprecated in dplyr 0.8.0.

We will re-write the code by using mutate_at() as shown in the code below.

col <- c("Country", "Continent")
globalPop <- read_xls("data/GlobalPopulation.xls",
                      sheet="Data") %>%
  mutate_at(col, as.factor) %>%
  mutate(Year = as.integer(Year))

Instead of using mutate_at(), across() can be used to derive the same outputs.

col <- c("Country", "Continent")
globalPop <- read_xls("data/GlobalPopulation.xls",
                      sheet="Data") %>%
  mutate(across(col, as.factor)) %>%
  mutate(Year = as.integer(Year))

3.3 Animated Data Visualisation: gganimate methods

gganimate extends the grammar of graphics as implemented by ggplot2 to include the description of animation. It does this by providing a range of new grammar classes that can be added to the plot object in order to customise how it should change with time.

• transition_*() defines how the data should be spread out and how it relates to itself across time.

• view_*() defines how the positional scales should change along the animation.

• shadow_*() defines how data from other points in time should be presented in the given point in time.

• enter_*()/exit_*() defines how new data should appear and how old data should disappear during the course of the animation.

• ease_aes() defines how different aesthetics should be eased during transitions.

3.3.1 Building a static population bubble plot

In the code below, the basic ggplot2 functions are used to create a static bubble plot.

ggplot(globalPop, aes(x = Old, y = Young, 
                      size = Population, 
                      colour = Country)) +
  geom_point(alpha = 0.7, 
             show.legend = FALSE) +
  scale_colour_manual(values = country_colors) +
  scale_size(range = c(2, 12)) +
  labs(title = 'Year: {frame_time}', 
       x = '% Aged', 
       y = '% Young') 

3.3.2 Building the animated bubble plot

In the code below:-

• transition_time() of gganimate is used to create transition through distinct states in time (i.e. Year).

• ease_aes() is used to control easing of aesthetics. The default is linear. Other methods are: quadratic, cubic, quartic, quintic, sine, circular, exponential, elastic, back, and bounce.

ggplot(globalPop, aes(x = Old, y = Young, 
                      size = Population, 
                      colour = Country)) +
  geom_point(alpha = 0.7, 
             show.legend = FALSE) +
  scale_colour_manual(values = country_colors) +
  scale_size(range = c(2, 12)) +
  labs(title = 'Year: {frame_time}', 
       x = '% Aged', 
       y = '% Young') +
  transition_time(Year) +       
  ease_aes('linear')          

3.4 Animated Data Visualization: plotly

In Plotly R package, both ggplotly() and plot_ly() support key frame animations through the frame argument/aesthetic. They also support an ids argument/aesthetic to ensure smooth transitions between objects with the same id (which helps facilitate object constancy).

3.4.1 Building an animated bubble plot: ggplot() method

Here we will learn how to create an animated bubble plot by using ggplotly() method.

The plot

The animated bubble plot above includes a play/pause button and a slider component to control the animation.

The code

gg <- ggplot(globalPop, 
       aes(x = Old, 
           y = Young, 
           size = Population, 
           colour = Country)) +
  geom_point(aes(size = Population,
                 frame = Year),
             alpha = 0.7, 
             show.legend = FALSE) +
  scale_colour_manual(values = country_colors) +
  scale_size(range = c(2, 12)) +
  labs(x = '% Aged', 
       y = '% Young')

ggplotly(gg)

Note

• Appropriate ggplot2 functions are used to create a static bubble plot. The output is then saved as an R object called gg.

• ggplotly() is then used to convert the R graphic object into an animated svg object.

Although show.legend = FALSE argument was used, the legend still appears on the plot. To overcome this problem, theme(legend.position='none') should be used as shown in the plot and code chunk below.

The plot

The code

gg <- ggplot(globalPop, 
       aes(x = Old, 
           y = Young, 
           size = Population, 
           colour = Country)) +
  geom_point(aes(size = Population,
                 frame = Year),
             alpha = 0.7) +
  scale_colour_manual(values = country_colors) +
  scale_size(range = c(2, 12)) +
  labs(x = '% Aged', 
       y = '% Young') + 
  theme(legend.position='none')

ggplotly(gg)

3.4.2 Building an animated bubble plot: plot_ly() method

In this section, we will learn how to create an animated bubble plot using the plot_ly() method.

The plot

The code

bp <- globalPop %>%
  plot_ly(x = ~Old, 
          y = ~Young, 
          size = ~Population, 
          color = ~Continent,
          sizes = c(2, 100),
          frame = ~Year, 
          text = ~Country, 
          hoverinfo = "text",
          type = 'scatter',
          mode = 'markers'
          ) %>%
  layout(showlegend = FALSE)
bp

3.5 References

• Getting Started

• Visit this link for a very interesting implementation of gganimate by your senior.

• Building an animation step-by-step with gganimate.

• Creating a composite gif with multiple gganimate panels